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I> ' We consider neutralino dark matter within the framework of SUSY GUTs 

with non-universal gaugino masses. In particular we focus on the case of 
, SU(5) with a SUSY breaking F-term in the 1, 24, 75 and 200 dimensional 

representations. We discuss the 24 case in some detail, and show that the 
■ bulk dark matter region cannot be accessed. We then go on to consider the 

admixture of the singlet SUSY breaking F-term with one of the 24, 75 or 
200 dimensional F-terms, and show that in these cases it becomes possible 
to access the bulk regions corresponding to low fine-tuned dark matter. Our 
■j-j . results are presented in the (Mi,M2) plane for fixed M3 and so are useful 

rN I for considering general GUT models, as well as more general non-universal 

& ' gaugino models. 
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1 Introduction 



Supersymmetry (SUSY) at the TeV scale remains an attractive possibility for new 
physics beyond the Standard Model. SUSY helps in the unification of couplings 
in Grand Unified Theories (GUTs), and provides a resolution of some aspects of 
the hierarchy problem. In addition the lightest SUSY particle (LSP) may be a 
neutralino consisting of a linear combination of Bino, Wino and neutral Higgsinos, 
providing a consistent WIMP dark matter candidate [T] . For example the minimal 
super symmetric standard model (MSSM) with conserved R-parity provides such 
an LSP with a mass of order the electroweak scale. Although general arguments 
suggest that such a particle should provide a good dark matter candidate [2J, the 
successful regions of parameter space allowed by WMAP and collider constraints 
are now tightly restricted [3"]-[T£]. 

Such a restricted parameter space has lead to recent claims that supersymmetry 
must be fine-tuned to fit the observed dark matter relic density [TUj . This is 
a serious concern for supersymmetry, especially as much of the motivation for 
supersymmetry arises from fine-tuning arguments in the form of its solution to the 
hierarchy problem. In previous work [20]- [22] we quantitatively studied the fine- 
tuning cost of the primary dark matter regions within the MSSM. It was found that 
the majority of dark matter regions did indeed require some degree of fine-tuning, 
and that this fine-tuning could be directly related to the mechanism responsible for 
the annihilation of SUSY matter in the early universe that defined each region. The 
one region that exhibited no fine-tuning at all was the 'bulk region' in which the 
dominant annihilation mechanism is via t-channel slepton exchange. This region 
can be accessed in models in which the gauginos have non-universal soft masses at 
the GUT scale @]. 

These results motivate a more careful study of models that give rise to non-universal 
gaugino masses. In our previous work such a region was accessed by allowing all the 
gaugino masses to vary independently. Such an approach is very unconstrained. 
We would expect the gaugino masses to arise from a deeper theory such as string 
constructions, as studied in [21], [22] or in GUT models [23]- [25]. Both approaches 
generally impose specific relations between the gaugino masses at the GUT scale. 
In this paper we shall discuss non-universal gaugino masses in a more general way 
than previously, allowing for different relative signs of gaugino masses, focusing on 
SU (5) GUTs as an example, although it is clear that similar effects can be achieved 
in other GUTs such as SO( 10) or Pati-Salam. We shall show how the bulk region 
may be readily accessed in such models providing that the SUSY breaking sector 
arises from a combination of an SU(5) singlet 1, together with an admixture of 
one of the 24, 75 or 200 representations of SU(5). We will also show that in all 
cases the fine-tuning required to access such a region remains small. 
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The rest of the paper is set out as follows. First we review our methodology 
in section [2J In section [3] we review the structure of gaugino non- universality in 
SU(5). In section H] we consider the specific case where all of the gaugino masses 
arise from a 24 of SU(5). In section [5] we generalise this to the case where the 
masses arise from an admixture of the singlet representation and one of the 24, 75 
or 200. In section [U] we present our conclusions. 

2 Methodology 

2.1 Codes 

The GUT structure of the theory is a structure that is imposed on the soft SUSY 
breaking masses at the GUT scale, incur ~ 2 x 10 16 GeV. To study the low energy 
phenomenology of such a model we need to run the mass spectrum down to the 
electroweak scale. To do this we use the RGE code SoftSusy [26J. This interfaces 
with the MSSM package within micrOMEGAs [27]. We use this to calculate the dark 
matter relic density QcDMh 2 , as well as BR{b — > sj) and 5a M . 

2.2 Experimental Bounds 

Not all choices of parameters are equal. After running the mass spectrum of the 
model point from the GUT scale to the electroweak scale we perform a number of 
checks. A point is ruled out if it: 

1. doesn't provide radiative electroweak symmetry breaking (REWSB). 

2. violates mass bounds on particles from the Tevatron and LEP2. 

3. results in a lightest supersymmetric particle (LSP) that is not the lightest 
neutralino. 

In the remaining parameter space we plot regions that fit BR(b — > 57) and Sa^ at 
la and 2a. 

2.2.1 5a M 

Present measurements of the value of the anomalous magnetic moment of the muon 
a M deviate from the theoretical calculation of the SM valud^. Taking the current 

1 There is a long running debate as to whether the calculation of the hadronic vacuum po- 
larisation in the Standard Model should be done with the e + e~ data, or the r. The weight of 
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experimental world average, and state of the art Standard Model value from [28] 
there is a discrepancy: 

K) e:cp - (a^sM = Sa„ = (2.95 ± 0.88) x 1(T 9 (2.1) 

which amounts to a 3. 4a deviation from the Standard Model value. 

2.2.2 BR{b -> si) 

The variation of BR(b — > sj) from the value predicted by the Standard Model 
is highly sensitive to SUSY contributions arising from charged Higgs-top loops 
and chargino-stop loops. To date no deviation from the Standard Model has been 
detected. We take the current world average from [2H] of the BELLE [3D], CLEO 
[31] and BaBar [32] experiments: 

BR(b -> S7) = (3.55 ± 0.26) x 10" 4 (2.2) 

2.2.3 QcDhih 2 

Evidence from the CMB and rotation curves of galaxies both point to a large 
amount of cold non-baryonic dark matter in the universe. The present measure- 
ments [33] place the dark matter density at: 

VLcDhih 2 = 0.106 ± 0.008 (2.3) 



For any point that lies within the 2a allowed region we calculate the fine-tuning 
and plot the resulting colour-coded point. 

2.3 Fine-tuning 



As in [zDJ we follow Ellis and Olive [M] in quantifying the fine-tuning price of 
fitting dark matter with the measure: 



A f 



d In (ttcDNfh 2 



d In (a) 



[2.4) 



where we take the total fine-tuning of a point to be equal to the largest individual 
tuning, A = max(A a ). 



evidence indicates the e + e data is more reliable and we use this in our work. 
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3 Gaugino Non-universality in SU(b) 



In the non-universal SU(5) model [16] , in addition to the singlet F-term SUSY 
breaking, the gauge kinetic function can also depend on a non-singlet chiral super- 
field $, whose auxiliary F-component acquires a large vacuum expectation value 
(vev). In general the gaugino masses come from the following dimension five term 
in the Lagrangian: 

< >ij 



L 



M 



Planck 



-XiXj 



(3.5) 



where are the U(l), SU{2) and SU(3) gaugino fields i.e. the bino B, the wino 
W and the gluino g respectively. Since the gauginos belong to the adjoint repre- 
sentation of SU(5), $ and can belong to any of the irreducible representations 
appearing in their symmetric product, i.e. 



(24 x 24) 



symm 



1 + 24 + 75 + 200 



(3.6) 



The minimal supergravity (mSUGRA) model assumes $ to be a singlet, which 
implies equal gaugino masses at the GUT scale. On the other hand if $ belongs 
to one of the non-singlet representations of SU(5), then these gaugino masses are 
unequal but related to one another via the representation invariants. Thus the 
three gaugino masses at the GUT scale in a given representation n are determined 
in terms of a single SUSY breaking mass parameter by 



Mi 



2,3 



^l,2,3 m l/2 



(3.7) 



where C\ 2 



.3 



1,1,1] 



Cf 



,2,3 



-3,2] 



Cf 



2.3 



-5,3,1) and OT 3 



'10, 2, 1). The resulting ratios of Mj's for each n are listed in Table HJ Of course in 



n 


M 3 


M 2 


Mi 


1 


1 


1 


1 


24 


1 


-3/2 


-1/2 


75 


1 


3 


-5 


200 


1 


2 


10 



Table 1: Relative values of the SU(3), SU(2) and U(l) gaugino masses at GUT 
scale for different representations n of the chiral superfield $. 



general the gauge kinetic function can involve several chiral superfields belonging 
to different representations of SU(5) which gives us the freedom to vary mass ra- 
tios continuously. In this, more general, case we can parameterise the GUT scale 
gaugino masses as: 

Mi, 2 ,3 = Cl Vi m n 1/2 (3.8) 
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where rriy 2 is the soft gaugino mass arising from the F-term vev in the represen- 
tation n. 

These non-universal gaugino mass models are known to be consistent with the 
observed universality of the gauge couplings at the GUT scale [23]- [2S], [3S] 

«3 = «2 = oil = «(— 1/25) (3.9) 

Since the gaugino masses evolve like the gauge couplings at one loop level of the 
renormalisation group equations (RGE), the three gaugino masses at the elec- 
troweak scale are proportional to the corresponding gauge couplings, i.e. 

(a 1 /a G )M 1 ~ (25/60) C?m? /2 
(a 2 /a G )M 2 ~ (25/30)C 2 n m? /2 
(a 3 /a G )M 3 ~ (25/9)C>? /2 (3.10) 

For simplicity we shall assume a universal SUSY breaking scalar mass m at the 
GUT scale. Then the corresponding scalar masses at the EW scale are given by 
the renormalisation group evolution formulae [36J . 



Mf = 
M 2 EW = 
Mi w = 



4 The 24 model 

We have previously seen [20] that a ratio Mi : M 2 : M 3 = 0.5 : 1 : 1 allows us 
to access the bulk region without violating LEP bounds. The bulk region in the 
CMSSM is usually ruled out because of a light Higgs. By allowing M 3 to be large 
we can avoid a light Higgs while allowing M\ to be light enough to give a light bino 
neutralino and light sleptons. This enhances neutralino decay via light t-channel 
slepton exchange and gives access to the bulk region. 

From Tabled] we observe that only the 24 model predicts a mass ratio M\ < M 3 . 
Therefore we shall explore the 24 model first. For the 24 model we have the input 
parameters: 

aE {m , my 2 , A , tan/3, sign(/i)}. 

where the masses are all set as in the CMSSM except for the gaugino masses which 
have the form: 

M x = -0.5 mf /2 
M 2 = -1.5 mf /2 
M 3 = m\% 

With this gaugino mass structure, the bino mass in the 24 for a given m\n is half 
of the bino mass in the CMSSM for the same mi/ 2 - The bino mass also affects the 
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running of the slepton masses such that lower Mi corresponds to a lower slepton 
mass. Therefore the 24 will have lower mass sleptons than the CMSSM for a given 
value of mo and mi/2- Light sleptons enhance the annihilation of neutralinos via 
t-channel slepton exchange (giving rise to a WMAP region known as the bulk 
region). Therefore we expect the bulk region to appear at larger ro^j than in the 
CMSSM and thus circumvent the Higgs mass bound. 



|— I 3 

CMSSM with A = 0, tan 3 = 10, sign(/i) = l. SU(5), n=24 with A = 0, tan 8 = 10, sign(//) = l. 




Figure 1: The parameter space for the CMSSM (top- left), the 24 model with sign(ju) +ve 
(top-right) and with sign(/x) — ve (bottom). Low rriQ is ruled out as the f becomes the 
LSP(light green). Low is ruled out as rrih < HlGeV. In the remaining parameter 
space, the only strip of allowed dark matter is a f — Xi coannihilation strip which shows 
comparable degrees of tuning in all plots. 

To study this effect, we look at the (mo, mi/2) plane in both the CMSSM and the 
24 in Fig. [TJ The CMSSM is shown in the top-left panel, the 24 with fj, positive in 
the top-right panel and the 24 with fi negative is shown in the bottom-left panel. 
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In the CMSSM scan we can see that low m is ruled out as the stau becomes 
lighter than the neutralino. Low m\n is ruled out as rrif L < HlGeV. The contours 
of 1 and 2a for Sa^ (green short and long dashed lines respectively) are plotted 
in the remaining parameter space, showing that the current measurement of Sa^ 
favours low mo and m^. Finally the region that satisfies WMAP is plotted as a 
multicoloured strip that runs alongside the light green region ruled out by a stau 
LSP. This WMAP strip is mostly red. This colour coding refers to a log measure 
of the fine-tuning and can be read off via the log-scale on the right hand side. The 
tuning of the f coannihilation strip agrees with our previous findings. 

In the second and third panels of Fig. [T] we once again display the (mo, mi/2) 
plane but this time using the 24 model's soft gaugino masses with \x positive and 
negative respectively. In both cases, low mo is ruled out by a stau LSP and low 
mi/2 is ruled out by a light Higgs. 

The Sa^ and BR(b — > 57) values are significantly different in the 24 model than 
in the CMSSM. Firstly neither 24 plot has a region that agrees with the current 
measured value of Sa^ (they both give 5a^ ± O(10 -10 )). Secondly BR{b — > 57) 
becomes an important constraint. For \i +ve, the model agrees with the measured 
value of BR{b — > sj) at la for large mi/ 2 (> 700 GeV) and agrees at 2a for low 
rax j2- With fi — ve, only the parameter space at m > 700 GeV fits BR{b — > 57) 
at 2a. Lower m exceeds this limit. 

Now consider the change in the dark matter strip. We expected to be able to 
access the bulk region in this model as we would have a lighter bino neutralino 
and lighter sleptons in the 24 model than in the CMSSM. This should move the 
bulk region to larger values of mi/ 2 and out from under the region ruled out by 
the LEP2 bound on the lightest Higgs boson. 

Contrary to our naive expectations, though the bulk region has moved to larger 
mi/2 in the 24 model, it remains ruled out. This is because the gaugino mass 
relations in the 24 also result in a lighter Higgs mass than the CMSSM, for the 
same m , m^. The only difference between the CMSSM and the 24 model is the 
magnitude and sign of the Mi and M 2 gaugino masses. Therefore the Higgs mass 
must be sensitive either to the sign difference between M 12 and M 3 or the larger 
value of M 2 . 

First consider the effect of the relative sign between M 12 and M 3 . In most RGEs 
the gaugino masses appear squared, however the trilinear RGEs have the form: 

dA t : 



dt 8vr 



2 



-1 ry 1 

tfAt + \Y b \ 2 A b + ( —g 2 3 M 3 + 3g 2 2 M 2 + — &M X 



(4.1i; 



If all Mi are positive, then the gauginos provide a large positive contribution to 
the RGE and so help to push the trilinear negative through the running. This in 
turn affects the running of the Higgs mass. In the 24 case, the sign of Mi 2 are 
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Figure 2: Here we show the running of At from the GUT scale value of At = to the 
weak scale for the point tuq = lOOGeV, m l / 2 = 350GeV, tan/3 = 10, Aq = 0. The 
running for the CMSSM is shown in blue, the running for the 24 model is shown in red. 



opposite to that of M 3 and so they reduce the contribution from the Gauginos 
and thus reduce the magnitude of the running, resulting in a small trilinear at 
the electroweak scale. Now we note that the contribution of M\ t2 are suppressed 
relative to that of M3 by a factor of gf, but this is partially compensated by the 
fact that |M 2 | > |M 3 | at the GUT scale. Therefore both the sign and magnitude of 
M 2 (GUT) are responsible for a substantial change in the running of the trilinears. 
This is shown in Fig. [21 

The change in the trilinear affects the running of m 2 H via the RGE: 



dm 2 H 
dt 



in" 



3\Y t \ 2 (m 2 Q:i + m 2 U3 + m 2 Hu + \A t \ 2 ) 
2 3 



-(3g 2 2 \M 2 



:9t \M 1 



(4.12) 



A smaller top trilinear results in a smaller running of the Higgs mass and a lighter 
Higgs. Therefore, as the 24 model results in a smaller value of A t at all energies 
below the GUT scale, it gives a smaller mass for the lightest Higgs than for the same 
model point in the CMSSM. This means that the LEP mass bounds for the lightest 
Higgs are more restrictive in the 24 model than in the CMSSM. Unfortunately, this 
results in the LEP Higgs bound ruling out the bulk region for all interesting regions 
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of parameter space of the 24 model. 



5 Two SU(b) Sectors 

We have seen that neither the CMSSM, corresponding to a singlet SUSY breaking 
sector, nor the 24 model is capable of accessing the bulk region of neutralino 
parameter space. Equally, as the 75 and 200 models have \M\\ > \Ms\, these 
sectors are even worse. In this section we therefore consider the next simplest 
possibility, namely that of two different SUSY breaking 577(5) representations 
acting together. Indeed, once one has accepted the existence of a single 24, 75 
or 200 dimensional SUSY breaking sector, it seems perfectly natural to allow the 
standard singlet SUSY breaking sector at the same time. In practice it may be 
difficult to avoid this scenario. 

Therefore we shall focus on the three simplest scenarios. We take the cases of a 
SUSY breaking sector consisting of: 

A (1 + 24) 
B (1 + 75) 
C (1 + 200) 

If we were to extend our model to allow three or four S77(5) representations con- 
tributing to SUSY breaking at once, we would be able to produce any pattern of 
non-universal gaugino masses. By constraining our model to two sectors we provide 
restrictions on the choice of gaugino masses which makes access to the bulk region 
non-trivial, and provides insight into what ingredients are required to achieve it. 



Mass 


A (1 + 24) 


B (1 + 75) 


C (1 + 200) 


Mi 
M 2 
M 3 


m\ /2 - 0.5 mf /2 
m\ /2 - 1.5 m\ A /2 

m l/2 + m l% 


m\ /2 - 5 ml% 

m \/2 + 3 m l% 
m \/2 + m l% 


m\ /2 + 10 m*° 
m \/2 + 2 m*f 2 
m \/2 + m T/°2 



Table 2: The gaugino mass relations for the different (1 + n) SUSY breaking 
scenarios. 

Within these models, we have different gaugino mass relations, shown in Table. [21 
By varying the soft gaugino masses Tn{? 2 , we describe three planes in the M lj2 ,3 
parameter space. 

Our aim is to access the bulk region. In [20] we found that the bulk region can 
be accessed in a model with non-universal gaugino masses for mo = 50 — 80 GeV. 
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Mi (GeV) Mi(GaV) 



1 - 


- 24 mixture 


wrong LSP 




>- 75 mixture 


■ no REWSB 




- 200 mixture I 


LEP mass limit 



Figure 3: The (Mi, M<i) plane with non-universal gaugino masses defined at the GUT 
scale. We take mo = 70 GeV, Aq = and tan/3 = 10 throughout vary M 3 : (a) M 3 = 
300 GeV, (b) M 3 = 400 GeV, (c) M 3 = 500 GeV, (d) M 3 = 600 GeV. For fixed M 3 , the 
allowed parameter space for each GUT mixture plotted as a line the (Mi, M2) parameter 
space. The WMAP allowed regions correspond to the elliptical regions in each quadrant, 
and are partially obscured by disallowed regions in panels (a) and (b). The BR(b — > S7) 
and 5a ^ regions are displayed as in Fig. [Hand discussed in the text. 

Therefore we fix m = 70 GeV, A = and tan/3 = 10. In Figs. E]^a)-(d) we plot 
the (Mi, M 2 ) plane for increasing values of M 3 , from 300 — 600 GeV. As Mi and 
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M 2 can in general be either positive or negative in (1 + n) scenarios, we allow Mi 
and M 2 to take positive and negative values. For a given M 3 , the gaugino mass 
relation of Table [2] constrain each of the (1 + n) scenarios to a line in the (Mi, M 2 ) 
plane. We plot these lines for each case. 

As each model has the singlet representation as a limit when m™/ 2 — > 0, all the 
lines converge at a point. At this point the model is precisely that of the CMSSM, 
and as such is ruled out for almost all M3 by a f LSP or the LEP bound on 
the lightest Higgs. The other end of each line corresponds to the opposite limit 
m\ /2 = 0, m™ /2 = M 3 . 

We also plot the BR(b — > sj) and Sa^ constraints. The only region that doesn't fit 
BR(b — > S7) within 2a is panel (a) at large M 2 . The values of 5a^ are insensitive to 
M3. In the quadrant with Mi and M2 +ve we have the largest SUSY contribution 
to 5a^, enabling the model to fit 5a^ at la. In the quadrant with Mi +ve, M 2 -ve, 
the model can fit 5a^ at la. For negative Mi we get a negative SUSY contribution, 
5a fj,. If we were to plot the parameter space with \i negative, 5a ^ would have the 
opposite sign and the model would fit the observed value of 5a M for negative Mi. 

Finally, we plot the dark matter regions with colours corresponding to their fine- 
tuning calculated with respect to the general non-universal gaugino model with 
parameters: a G {mo, Mi, M 2 , M3, Aq, tan/3}. This allows us to easily pick 
out the bulk region as it is 'supernatural' with A n < 1 and is therefore plotted in 
yellow. We use this to pick out the points at which each (1 + n) representation 
provides access to the bulk region. We take these points and calculate the dark 
matter fine-tuning with respect to the (1 + n) model in question. 

First consider the 1 + 24 model. In Figs. E^a), (b) the model does not access the 
bulk region. This fits with our results of section H] as low mi/2 is ruled out by a 
light Higgs in the 24 scenario. In Figs. EJ^c), (d), we can access the bulk region with 
a mixture that is primarily 24. We show the corresponding fine-tuning for both 
points in Table El Note that for both points rriy 2 > rn\, 2 , so the gaugino masses 
arise predominantly from the 24. 

Next consider the 1 + 75 model. This model lies along the blue short dashed line. 
The 75 limit is not shown. This is because in the pure 75 scenario Mi = — 5M 3 . 
Therefore the 75 limit lies outside the range plotted for all M3 that we consider. 
In such a limit, as studied in p3], [16], the lightest neutralino is predominantly 
higgsino. As discussed earlier we cannot access the bulk region in such a limit. 
This limit lies off the plots and we do not consider it further here. 

In the 75, Mi is negative. This results in two scenarios in which Mi < M 3 . For a 
small Tn[ 5 , 2 , the negative contribution results in a small, positive, Mi. For a slightly 
larger 7Ttw 2 , we get a small, negative Mi. This is shown in the plots and is the 
reason that the 1 + 75 accesses the bulk region twice for most values of M 3 , once 
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Parameter 


Al 


A2 




value 




value 




m 


70 


1.43 


70 


0.96 


m i/2 


33.3 


0.026 


100 


0.39 


1/2 


466.7 


0.075 


500 


1.02 


A Q 














tan [3 


10 


0.37 


10 


0.21 


Max 




1.43 




0.96 


Mi 


-200 


0.19 


-150 


0.59 


M 2 


-666.7 


0.21 


-650 


0.38 


M 3 


500 


0.075 


600 


0.0088 



Table 3: The fine-tuning for points Al and A2 that lie within the bulk region for the 
(1 + 24) model. For both points Triy 2 > rn\i 2 , so the gaugino masses arise predominantly 
from the 24. In the lower section of the table we give the corresponding GUT scale Mj for 
each point. As the tunings plotted in Fig. [3] are calculated with respect to the parameter 
set a G {mo, Mi, Mi, Ms, Aq, tan/3}, we give the relevant tunings with respect to the 
individual M, for comparison. 

for each sign of Mi. We study the 7 resulting points in the bulk regions in Table HI 
Note that for all points my 2 < rn{/ 2 , so the gaugino masses arise predominantly 
from the singlet. 

Finally consider the case of the 1 + 200 model. The lines corresponding to this 
model are plotted in red with long dashes. As in the 1 + 75 case, in the 200 limit 
the lightest neutralino is higgsino and we cannot access the bulk region. This limit 
lies off the plots and we do not consider it further here. 

As the 200 has all gaugino masses positive, and large M\, we cannot access the 
bulk region in the 200 limit. However by combining with the singlet we can get 
I Mi I < |M 3 | by taking a small, negative mfj 2 . This allows such a model to access 
the bulk region for positive and negative small Mi. We study the resulting 6 points 
in the bulk region in Table [51 In all points l^^l < l m i/2l so the gaugino masses 
arise predominantly from the 1. 

The hierarchy of the weak scale SUSY spectrum is fairly stable for all the points 
shown in Fig [3j Table [6] lists the neutralino, chargino and sfermion masses along 
with Mi, M 2 and the Higgsino mass parameter \x for the point B5 as an example. 
In contrast to the CMSSM the bino is lighter than the wino by a factor of 6. 
Correspondingly the right and left slepton masses are split by a large factor. The 
small value of mo also ensures that the right handed sleptons are considerably 
lighter than the wino. Hence a large fraction of wino decay is predicted to proceed 
via fi, resulting in one or more tau leptons in the final state in addition to the 
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Parameter 


Bl 


B2 


B3 


B4 




value 




value 




value 




value 




m 


70 


0.91 


70 


1.18 


70 


0.86 


70 


1.0 


m \/2 


217 


0.78 


300 


0.64 


363 


1.4 


387 


1.1 


7^ 

771 1/9 
1/2 


83.3 


1.4 


100 


0.91 


36.7 


0.67 


113 


1.5 




























tan /3 


10 


0.13 


10 


0.29 


10 


0.14 


10 


0.32 


Max 




1.4 




0.91 




1.4 




1.5 


Mi 


-200 


0.66 


-200 


0.38 


180 


0.67 


-180 


0.51 


M 2 


467 


0.086 


600 


0.032 


473 


0.096 


727 


0.075 


M 3 


300 


0.13 


400 


0.071 


400 


0.061 


500 


0.047 


Parameter 


B5 


B6 


B7 








value 


A n 


value 


A u 


value 


A u 






m 


70 


0.75 


70 


0.95 


70 


0.84 






m i/2 


450 


1.8 


475 


1.7 


530 


2.0 






7^ 

m 1/9 
1/2 


50 


0.99 


125 


2.4 


70 


1.2 






^0 
























tan /3 


10 


0.15 


10 


0.32 


10 


0.22 






Max 




1.8 




2.4 




2.0 






Mi 


200 


0.80 


-150 


0.55 


180 


0.64 






M 2 


600 


0.038 


850 


0.082 


740 


0.031 






M 3 


500 


0.014 


600 


0.16 


600 


0.12 







Table 4: The fine-tuning for points Bl-7 that lie within the bulk region for the (1 + 75) 
model. For all points niy 2 < iriy 2 , so the gaugino masses arise predominantly from the 
singlet. In the lower section of the table we give the corresponding GUT scale Mj for 
each point. As the tunings plotted in Fig.[3]are calculated with respect to the parameter 
set a € {mo, Mi, M 2 , M3, Aq, tan/3}, we give the relevant tunings with respect to the 
individual Mj for comparison. 

missing- Et- Though the light selectron and smuon have negligible left-handed 
components, and so cannot take part in the wino decay, the heavier selectron and 
smuon are still lighter than the wino in all points we consider. A wino decay via 
a left-handed selectron/smuon would give a distinctive signal in the form of hard 
electron (s)/muon(s) in addition to the missing- Et- Thus one expects a distinctive 
SUSY signal from squark/gluino cascade decays at LHC containing hard isolated 
leptons in addition to the missing- Et and jets. 
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Parameter 


CI 


C2 


C3 




value 




value 




value 




m 


70 


1.6 


70 


0.89 


70 


1.1 


m i/2 


467 


0.11 


424 


1.7 


576 


1.4 


^200 

III-, In 

L/Z 


-66.7 


0.40 


-24.4 


0.93 


-75.6 


2.2 


A 




















tan (3 


10 


0.79 


10 


0.25 


10 


0.54 


Max 




1.6 




1.7 




2.2 


M x 


-200 


0.19 


180 


0.67 


-180 


0.56 


M 2 


333 


0.83 


376 


0.31 


424 


0.59 


M 3 


400 


0.75 


400 


0.22 


500 


0.39 


Parameter 


C4 


C5 


C6 




value 


A u 


value 


A n 


value 


A n 


m 


70 


0.78 


70 


0.97 


70 


0.86 


m i/2 


533 


2.3 


683 


2.3 


647 


2.5 


^200 

llb-i In 

1/2 


-33.3 


1.3 


-83.3 


3.1 


-46.7 


1.7 






















tan /3 


10 


0.23 


10 


0.43 


10 


0.25 


Max 




2.3 




3.1 




2.5 


M x 


200 


0.80 


-150 


0.59 


180 


0.63 


M 2 


467 


0.25 


517 


0.49 


553 


0.22 


M 3 


500 


0.13 


600 


0.20 


600 


0.047 



Table 5: The fine-tuning for points Cl-6 that lie within the bulk region for the (1 + 200) 
model. For all points Ijto^I ^ l 771 !^!' so ^ e g au gi no masses arise predominantly from 
the 1. We also give the corresponding GUT scale Mj for each point. As the tunings in 
Fig. [3] are calculated with respect to the parameters a £ {mo, M 2 , M3, Aq, tan/3}, 
we give the tunings with respect to Mj for comparison. 

6 Conclusions 

In previous work we found that a model with non-universal gaugino masses could 
access the bulk region in which t-channel slepton exchange alone could account for 
the observed dark matter relic density. The bulk region is an attractive prospect 
as it allows SUSY to account for the observed dark matter relic density without 
any appreciable fine-tuning. However, a model with entirely free gaugino masses 
is very unconstrained. Such non-universality must arise from a deeper structure 
and such structures should impose restrictions on the precise form of the gaugino 
masses at the GUT scale. 

In this paper we have considered neutralino dark matter within the framework of 
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Particle 


Mass (GeV) 


Xi ( bin o) 


78.1 


X2 (wmo) 


457 


X 3 (niggsino) 


614 


X4 (niggsino) 


636 


Xi (wmo) 


461 


Xj (niggsino) 


635 




81 


71 JEW 

Mf 


470 


/1/f-BVK 
M 3 


1 i on 

ilzO 




Oil 


9 


1150 


fi 


104 




399 




115 


eL,fiL 


399 




793 


k 


1025 




980 


k 


1000 


Ql,2,R 


~ 1005 


<2l,2,L 


~ 1070 



Table 6: The SUSY mass spectrum of point B5 from Fig. This spectrum is charac- 
teristic of all bulk region points we have studied. We display the hierarchy and flavour 
of the neutralino and chargino sectors. We also display the values of the neutralino mass 
parameters for completeness. For the squarks we take a typical squark mass rather than 
list the full squark spectrum. The exceptions are the 3rd family squarks that we list 
separately. Finally, the sneutrinos are degenerate with e,fitL- 
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SUSY GUTs with non-universal gaugino masses. We have taken the specific case 
of an 577(5) GUT model where the gaugino masses arise from different irreducible 
representations of the symmetric product of the adjoint representations. In partic- 
ular we focused on the case of SU (5) with a SUSY breaking F-term in the 1, 24, 75 
and 200 dimensional representations. We discussed the 24 case in some detail, and 
showed that the bulk dark matter region cannot be accessed in this case. In general 
if we just take the simplest case in which the gaugino masses arise from only one 
representation, we find that as far as achieving the bulk region is concerned, there 
is no advantage over the CMSSM. This is in part due to the surprising result that 
the sign and magnitude of M2 with respect to M3 has an important effect on the 
lightest Higgs mass through its effect on the top trilinear. 

We then went on to consider the case of the singlet SUSY breaking F-term com- 
bined with an admixture of one of the 24, 75 or 200 dimensional F-terms. Such 
a scenario is natural once we allow the higher dimensional representations in our 
theory. In all these cases we showed that it becomes possible to access the bulk 
regions corresponding to low fine-tuned dark matter. In addition, the degree of 
fine-tuning required to access the bulk region remains small in the GUT models. 
Therefore we conclude that such models can access the bulk region and naturally 
account for the observed dark matter relic density. 

Finally we note that the results in Fig. [3] are presented in the (Mi,M 2 ) plane for 
fixed M 3 and so are useful for considering general GUT models, as well as more 
general non-universal gaugino models. The hierarchy of weak scale SUSY spectrum 
is fairly stable for all the points shown in Fig. [3j Both the right and left sleptons are 
lighter than the wino, implying a large leptonic BR of wino decay. This promises 
a distinctive SUSY signal from squark/gluino cascade decays at LHC in the form 
of hard isolated leptons in addition to the missing- E T and jets. 
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